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Abstract 

We show that a four-dimensional equation of state for a cosmological constant term arises from 
a perfect fluid in the bulk in the context of a gravity model where the scalar curvature is non- 
minimally coupled to the perfect fluid Lagrangian density. The four-dimensional theory is fully 
determined from the induced equations on the brane, subject to the boundary conditions derived 
across the brane. 
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I. INTRODUCTION 



Braneworld scenarios are an interesting development in what concerns gravity models 
and their cosmological implications Most often these scenarios assume, given observa- 
tional constraints as well as theoretical assumptions, that the bulk space is empty except 
for a cosmological constant. However, more recently the implications of having vector and 
scalar fields in the bulk were studied in connection with Lorentz symmetry j^] and gauge 
symmetry breaking [sl]. 

In this paper, we consider the presence of a perfect fluid in the bulk in the context of 
a five-dimensional braneword model where the scalar curvature couples non-minimally to 
the Lagrangian density of the perfect fluid. In (3+1) dimensions this class of models with 
Lagrangian density of the form {4 1 

/: = ah{R) + {l + XMR))CM , (1) 

where fi{R) and f2{R) are generic functions of the scalar curvature, was shown to exhibit 
interesting features which allow one to address problems such as the rotation curves of 
galaxies without the need of dark matter (see Ref. and references therein) and the 
Pioneer anomaly (see Ref. 5|, l6|] and references therein). The stability of these models 
has been examined in Ref. j?]. Other studies on their implications included their impact 
on stellar equilibrium and the analysis of their corresponding PPN parameters, which were 
studied in Refs. , respectively. 

Recently there has also been interest in the conformal equivalence between f{R) theories 
and Einstein gravity non-minimally coupled to a scalar field in the context of braneworlds 
jlol | . The expected increasingly higher order of the discontinuity of the geometric quantities 
across the brane with the increasing power in R of f{R) is solved by enforcing continuity of 
the metric to correspondingly higher-order derivatives. Here, however, we will not impose 
further continuity conditions on the intervening fields, allowing for the discontinuity of the 
second derivative of the metric across the brane and orthogonal to its surface, despite also 
obtaining an increase in the power of R. 

Crucial in the setting of our problem is a suitable implementation of the Israel matching 
conditions in the presence of bulk fields in order to extract the boundary conditions, both 
for gravity and the matter fields, which the induced equations of motion on the brane must 



satisfy. The method to be employed here was first introduced in Ref. Uj and further 
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developed in Refs. [2|, [3']. For completeness and clarity, the more involved technical details 
of our method are presented in the Appendix. As we shall see, and rather remarkably, the 
projection of the bulk perfect fluid induces on the brane a new cosmological constant term. 
This new source for a brane cosmological constant opens quite interesting perspectives 
for inflation at the early universe and for acceleration at the late time expansion of the 
universe. This result suggests that a perfect fluid in the bulk space may have a bearing on 
the cosmological constant problem on the brane. 

This paper is organized as follows. In section [III we present our model and work out a 
suitable Lagrangian density for a perfect fluid. This development extends the approach of 
Hawking and Ellis to the bulk space. In section UTTl we work out the matching conditions 
across the brane and derive the equations of motion therein induced. A derivation of the 
Gauss-Codacci relations is also presented in the Appendix for completion. Section [IV] 
contains our results and section |V] our conclusions. 

II. A MODIFIED GRAVITY MODEL IN THE BULK 
A. The Einstein Equations 



We consider the particular case of the action discussed in Ref. j^. We set fi{R) = 
f2{R) = R and introduce a cosmological constant as follows 

S = j (fx^g [Ml,^,)R + (1 + XR)£.M - 2£a] . (2) 

Here Mp(5) is the five- dimensional Planck mass. Cm and Ca are respectively the matter 
and the cosmological constant Lagrangian densities. 
We define the stress-energy tensor as usual 

For convenience, we define also the vacuum energy tensor as 

assumed to be of the form 

^^.u = Ai5)9f.u + a5iN) {g^, - N,N^) , (5) 
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so as to include both the bulk vacuum value A(5) and the brane tension a. Here A*"^ are the 
components of the unit five- vector orthogonal to the brane = N^e^. Thus, in Gaussian 
coordinates, the cosmological constant tensor is given by: 



A 



(4) 



g(A(5)+CT) 



\ 



A 



(6) 



(5) / 



The five- dimensional Einstein equation is obtained by varying the action with respect to 
the metric, finding that 



1 



^p(5)G'm. - o (1 + T^.u + A^, - A (V^V, - g^^U - XR^,) C 



AI 







(7) 



B. A Perfect Fluid in the Bulk 



Since the Einstein equation in Eq. ([7]) contains terms in Cm, we must construct a 
Lagrangian density associated with a perfect fluid so that by Eq. ([3]) it yields 



T^,u = (p + p)u^Uy + pg^u 



(8) 



where p is the energy density, p the pressure and the unit five-velocity of the fiuid 
(tangent to the fiow lines and thus time-like, u^m^ = —1). For this purpose, we follow 
closely the procedure described in Ref. [l2|. 

Let the perfect fiuid Lagrangian density Cm be given by 



C 



M 



-r(l + e) 



(9) 



where r is an auxiliary variable and e is the internal energy of the fiuid as well as a 
function of r. Assuming that the fiuid current vector = tu^ is conserved, V^j'^ = 0, 
then S{y/—gj^) = when the metric is varied. Then, from 



it follows that 



r 



1 



-ffg^u = - iV^f) iV^f) g,u 



(10) 



2t5t 



= {j,Ju-fjpg,u)Sg^'' , (11) 
and consequently that the variation of r with respect to the metric is given by 

^'^ = 1 ('^9,.u + ruf,u^) Sg^"" . (12) 

Using the definition of the stress-energy tensor, Eq. ([3]), we find that 

.5Cm 



T^u — ^Mgfiu ~ 2- 



r(l + e) + T^ 



2 de 



dr 

ip + p)ui,u^+pg^u , (13) 
where we made the following identifications 

p = T{l + e), P = ^^%- (14) 

We have thus obtained the stress-energy tensor for a perfect fiuid from the Lagrangian 
density Cm = —p and from the continuity equation V^(rM'^) = 0. For an alternative 
formu 
Refs. 

To obtain the equation of motion for the perfect fiuid, one could compute the diver- 
gence of the Einstein equation, Eq. ([7]). Alternatively, we will proceed directly from the 
Lagrangian density. For this purpose, we consider the action of the Lie derivatives on the 
fiuid fiow lines. Let 7 : [a, 6] x A/" ^ P C be a congruence of fiow lines, one through 
each point of Ai, where [a, b] is the interval of the parameter ascribed to the fiow lines, Af is 
some four-dimensional manifold and P is a small region of the five-dimensional spacetime 
manifold Ai. The tangent vector to the fiow lines is given by U = {d/dt)^, with t G [a, b], 
which once normalized we identify with the fiuid velocity 



ation where the Lagrangian density is identified with the pressure see, for instance, 

3,Q. 



u'' = , = iV7i ■ (15) 

V-9apU'^UP \U\ 

The action S is required to be stationary when the fiow lines are varied. Variations with 
respect to the fiow lines, which we here represent by A, amount to variations along the 
corresponding tangent vectors. Considering 7(r, [a,6],M'^), where r is the parameter that 
selects different congruences of fiow lines, then AU = £vU, where V = {d/dr)^ and £v 
is the Lie derivative along V. Since A(M^|f/|) = C^U, then 

AM^ = ^(/:vf/^-w^A|f/|) . (16) 



Moreover, 



and 



and consequently 



A\U\ = -^(?„^f/"Af/^ = -g^fsu'^C^U^ 



(17) 



From the conservation of the fluid current vector j^.^ = 0, it follows that A{j^. 
(Aj^);^ which, using Eq. (JT9l) and integrating along the flow lines, yields 

At = {TVf^)-p + tV^.^u'^u'' . 



(19) 
= 

(20) 



Thus, in the Lagrangian density, r varies so that the associated current vector is con- 
served. 

Finally, the condition for the stationarity of the metric yields 
ds r . , (dC 



j y/^d^x + \R) 
r^d^x + \R) 



1 + 



d{re) 
dr 
. , djre) 
dr 



At 



[V, {rVn + r (V^V;) u^u"] | = . (21) 



Integrating by parts and using the Stokes theorem to discard the surface terms, we find 
that 



/-gd°x |rV^ 

+ V;3 



[l + XR) 



l + XR) (^1 + ^U„ ..P 



dr ' 



(22) 



which holds for any vector V. Therefore, the expression within curly brackets must vanish 



XVpR 



+ (l + XR) 



d{Te) 
dr 

_ ^ilS 
dr 



+ (1 + XR) 



d{ 



re 



dr 



T 



U^{VpU^')T = 



Using the equations in Eq. f|T4l) and noticing that 



T 



'djrey 
dr 



, de 
dr 



(23) 
(24) 

(25) 



we obtain the equation of motion for the perfect fluid in the bulk 



[A (p + p) VpR + (1 + \R)Vpp\ {g^^ + u^u^") + (1 + \R){p + p)u^V pu^ = . (26) 



For p + p 7^ 0, this equation can be rewritten as 



ds ds ' '""^ 
where can be regarded as an exterior force given by 



(27) 



1 + XR 



P + P 



(2^ 



Setting A = 0, one recovers the known equation of motion for a perfect fluid in General 
Relativity. Moreover, Eq. (1251) agrees with the result obtained in Ref. ^\ by taking the 
divergence of the Einstein's field equation. This indicates that the assumed conservation 
of the fluid current vector j'^ = tu'^ is a consistent description of our physical system. 

The continuity equation, which follows from the conservation of j^, provides the last 
equation and ensures that our problem is well defined. From Eq. (1141) we find that 

1 



P = re,f, — , P-^{P + p) = P,i.P ■ 
T,, 1 + e 



(29) 



The continuity equation 



is thus equivalent to 



P,f.u'' + {p + p)u^'.^ = . 



(30) 



(31) 



For the ideal fluid, the gravitational field enters only in Eq. ( 1261) . Gravity does not appear 
in the velocity equation, Eq. (!3Ti) . as the velocity is measured relative to freely moving 
observers [l5|. 



III. THE INDUCED EQUATIONS ON THE BRANE 

In this section, we derive the equations of motion induced on the brane. First, we rewrite 
the components of the equations of motion derived in the previous section in Gaussian 
normal coordinates. In our notation, the directions parallel to the brane are denoted by 
A,B,..., while the normal direction is denoted by N so that the brane is localized at 

7 



= 0. Using the results derived in Appendix |Xl we obtain the relevant components of the 
Einstein equation: 

M|(5) (^)G'ab - ^ (1 + a (^)i?) [(p + P)UAUB + PQAb] + ^AB 

+ A [VaVb + KabVn - 9AB (□ + KVn + Vat Vat) — ^^^Rab] P = , (32) 

+ A [VaVn - K^Vb - '^'^Ran] P = 0, (33) 

M|(5) '^'^Gnn - ^ (1 + ^ ^'^^) [(P + pX + P9nn] + A(5) 
+ A [VatVtvP - (□ + KVn + VtvV^) - ^''^Rnn] P = . (34) 

Note that we will only keep the index indicating the corresponding dimension on the five- 
dimensional terms and whenever confusion is propitiated. 

We then proceed to derive the matching conditions, which follow from the presence of 
the brane dividing into two regions the bulk spacetime and from the symmetries of the 
bulk fields across the two regions. Here we regard the brane as a Z2-symmetric surface of 
infinitesimal thickness 26 and thus separating the bulk into two mirroring regions about 
= 0. The symmetry about the brane establishes how bulk quantities relate on the 
two sides of the brane. Hence, vector components parallel to the brane are even in A^, 
whereas normal components are odd. For tensor quantities this generalizes by considering 
each additional normal component to reverse the parity of the component with one less 
normal component. Accordingly, it follows that ua{N = —6) = ua{N = +6), whereas 
u,^{N = -6) = -un{N = +6). Likewise, qab^N = -6) = qab^N = +6) and Kab{N = 
—6) = —Kab{N = +6). Consequently, there will be quantities that are discontinuous across 
the brane and whose derivatives in A^ generate singular distributions at the position of the 
brane. Integration of these contributions in the coordinate normal to the brane allows to 
relate the induced geometry of the brane with the induced stress-energy therein localized. 
Hence, by extracting the singular contributions from the projected bulk equations and 
establishing the matching conditions, we obtain the equations of motion induced on the 
brane. 

We equate the (AB) component of the Einstein equation using the Gauss-Codacci con- 
ditions in Eqs. flA201lA22l) as well as Eq. flA13l) . Integrating along the A^ direction across 
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the position of the brane, we find that 



h5 

- agAB = lim / rfiVV^<' M'L^. {-Kab + QabK) 



5^0 



-5 



+\K [{p + p)uaUb + PQab] - A {qab^n - Kab) P | , (35) 

where upon integration by parts non-singular terms arise which vanish upon integration and 
which contribute to the effective equation of motion. For that we assume the energy density, 
the pressure and the fluid velocity to be continuous. This implies that only higher than 
second order derivatives in the direction can be singular on the brane and consequently 
survive over the inflnitesimal integration. From these considerations there follows the Israel 
matching condition 

(Mp(5) - Xp) {-Kab + QabK) + \K{p + p) {uaUb + Qab) - X^Qab^nP = -^(^Qab (36) 
which, upon taking the trace, yields 

(M|(5) - Xp) K{-1 +d) + XK{p + p) {ul + d)- XdVNp = -^CT . (37) 
Another useful result is 

(M|(5) - Xp) Kab = K^ [qab {mI^^^ - Xp) + X{p + p) [duAUB - QabuI)] . (38) 
Substituting Eq. (l36l) back in the {AB) component of the Einstein equation, we obtain 

+ 9AB^{5) 

[{p + p) {uaUb + Qab)] 
p = . (39) 



(M|(5) - Xp) 
1 



Gab - KKab + \gAB [K^ + KcdK'''') 



+ 



-{l + X{R-K^ - KcdK'"']) + XKVn 



^gAB + X (VaVb - qab^'c) 



From the {AN) component of the tensor equation we notice that 

Gan = Kab\b - K\A = -Vb J dNGAB = -V aTab , (40) 

where Tab stands for the stress-energy tensor induced on the brane as given by the Israel 
matching condition in Eq. fl36|) . If we impose conservation of energy on the brane, it follows 
that Gan = 0, which implies the condition 

Va [K{p + p){uaUb + gAs)] - [gAB^N - {K AB - gABK)] V Ap = • (41) 
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Furthermore, equating the [NN) component of the Einstein equation and integrating 
along the N direction, we find that 

h<5 



cINVn {pK + K[{p + p)u\ +p\] = -2K{p + p) 

-5 



C 



(42) 



Substituting back in the {NN) component of the Einstein equation, we find that 



(M|(5) - \p) - {-R + KcdK^^) + A(5) 



+ 



{1 + X{R-K^- KcdK^""]) + XKVn 



[{P + P) [nl + l)] 



- A (V^ + KVn) 



p = 0. 



(43) 



Moreover, substituting Vivp from the Israel matching condition for a time- like fiuid velocity 
normalized so that -u^ + = — 1, it follows that 



(M|,(5) - Xp) 



-R + K' 



d 



1 



+ 
+ 



{l + X{R-K^~ KcdK^''}) + XKV 



N 



+ A(5) 

[ip+pyc] 



p-^XK' (p + p) {ul + d)-^Ka = 0. 



(44) 



We treat analogously the equations of the perfect fiuid in the bulk. The equation of 
motion, Eq. (!26|l . can be combined with the continuity equation, Eq. (!3T|) . to yield 



V, [(1 + A (p + p) [g^, + u^u,)] - (1 + A (^)i?) g^,V,p = 0. 



(45) 



Substituting the expression for ^^^R in Eq. flAlOp . we integrate both the parallel and the 
orthogonal components along the direction to obtain the corresponding boundary con- 
ditions on the brane. From the parallel component we find that 

"+-5 



= dNVNf^-'2XVA[Kip + p)igAB + UAUB)] 

Kp + p) 



+XK {p + p)uaUn 



9ab\1 + -A+ . 

df d(Mf^^-Xp 



[duAUB - gABUc) 



+ [1 + X{R~K^- KcdK'''' - 2K^n)] (p + p)umUb + 2XKgAB'^AP^ ■ (46) 

Here we encounter a third derivative along the direction of the induced metric gAs- For 
a continuous metric, the first derivate can be discontinuous, the second derivative can be a 
delta-like singularity and consequently the third derivative can be a double-peaked delta. 
When we integrate in N we are left with a term in K at, which is proportional to the second 
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derivative of the metric and thus potentially singular, evaluated at the end points along the 
normal direction which define the thickness of the brane. Due to the Z2-symmetry, however, 
whereas the delta singularity is even about the brane, with K{N = —6) = —K{N = 
+6) and thus J dN'W nK = 2K, the double-peaked delta is odd, with K^fyr{N = —6) = 
+K^]\r{N = +6) and thus / dNV nK = 0. Consequently, when integrated along A^, only 
odd order derivatives of the metric along jump across the brane and thereby relating with 
the singular matter distribution at the location of the brane, while even order derivatives 
cancel at the end points. However, the term in question contains also the factor uaUn which 
is odd about the brane, thus causing the integral to survive and yield / dNV n (K^j^uaUn) = 
2K j^uaUn But N = +6. On the other hand, since the boundary condition in Eq. P2l) imposes 
that either ua = ot p + p = or K = 0, then regardless the case this term vanishes on 
the brane. Then, the parallel component of the boundary condition becomes 



2AVa [Kip + p) {qab + uaUb)] + 2\KgABVAp 



+ {p + p)ubUn 



. (47) 



Substituting the boundary condition in Eq. (H7|) back in the parallel projection of Eq. (H^ . 
and using also the energy conservation condition in Eq. (14T|) . we find for the induced 
equation for the fiuid on the brane 

[(1 + \{R~2K^- KcdK'''' + 2i^Viv}) (p + p) {qab + uaUb)] 
+ 2\K{p + p) [uaUb {Kab\c - Kac\b) + uaUn {KcdKdc9ab + KadKdb)] 

- XK^Va Hp + p) {uaub + Qab)] 

- XK ^AB ^1 + + djif—rx-p (^"^"^ - ^-^^"^) 

- [2\K {Kab - 9abK) + QAB (1 + \{R-K^- KcdK''''])] VaP = . (4J 



Vat [(p + p)uaUn\ 



The orthogonal component yields a trivial matching condition because of the continuity 
conditions across the brane of the quantities involved. We conclude that this component 
will only be important for the propagation of the fiuid off the brane and across the bulk. 
The propagation on the brane is solely described by Eq. (HHl) . which contains already the 
continuity condition in Eq. flHTl) . with the possibility of propagation off into the bulk being 
constrained by the conservation condition in Eq. (HTl) . 

These are the induced equations on the brane and can be solved the following way. The 
effective equations of motion in Eqs. (|39|1 and (j48l) consist of a coupled system which must 
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be solved together for the induced metric and for Viv [(p + p) {uaUb + Qab)] ■ With these 
results, we can then solve the stress-energy conservation condition in Eq. ( 14T1) . derived from 
the (^A^) component of the Einstein equation, and the (NN) component of the Einstein 
equation in Eq. (jH]) together for p and the extrinsic curvature, constrained by the matching 
conditions in Eqs. (!36|l and (H7|) which then allow to find the functional form of p in terms 
of p. 

From these equations we observe that the coupling of the curvature to the matter La- 
grangian density yields a contribution to the effective Newton's constant on the brane. 
Moreover, it is only through the non-minimal coupling that matter in the bulk interacts 
with that localized on the brane, in this case the tension a only. Notice that if A = 0, 
i.e. in the absence of a non-minimal coupling of the curvature to the matter Lagrangian 
density, Eqs. ( I47l) and ( l48l) read respectively 

p + p = (49) 

and 

Va(p + p) {qab + uaUb) - Qab^aP = . (50) 

This means that the presence of a perfect fluid in the bulk space induces on the brane an 
equation of state characteristic of a cosmological constant, without, however, the quantities 
p and p characterizing the fluid being necessarily constant. 

IV. RESULTS 

In this section we proceed to study the set of equations derived above for the three 
particular cases encompassed by the boundary condition derived from the (NN) component 
of the Einstein equation, Eq. (H2|) . 

A. Reflecting Boundary Condition: ua = 0, Vatutv = 

One of the cases that we can consider is that of the brane consisting of a reflecting surface 
for the incoming fluid flux from the bulk. This translates into setting Dirichlet boundary 
conditions to the parallel component of the fluid velocity ua = and Neumann boundary 
conditions to the normal component VatMat = at the position of the brane. Moreover, 
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the Z2-symmetry implies that VatMa = 0, whereas boost and translation invariance on the 
brane implies that V bUa = 0. 

The equations derived in the previous section become as follows. The (AB) and the 
(NN) components of the Einstein equation are given respectively by 



(M|(5) - \p) 
1 



Gab - KKab + ^9ab {KcdKcd + K^) 



+ 9AB^{5) 



+ 



{1 + X{R~K^- KcdKcd}) + XKVn 



(P + p)9ab 



Qab + A (VaVb - gAB^c) 



p = 



(51) 



(^P(5) - Ap) - 



-R + K' [--1] - KcdKcd 



+ A 



(5) 



+ 



i-AV?, 



p - \K^ (p + p)- -Ka = ; 



(52) 



the parallel component of the induced fluid equation and the condition of stress-energy 
conservation on the brane respectively by 

[(1 + A {i? - 2K^ - KcdK'''' + 2KVn}) (p + p)9ab] - XK^qab^aIp + p) 
- [2\K {Kab - 9abK) +gAB{l + X{R-K'^ KcdK''''})] VaP = 0, (53) 



^id-l) (M|(5) - Ap) VbK - A (^Kab - \9abK^ V aP = 



(54) 



We present also the matching conditions from both the {AB) component of the Einstein 
equation, Eq. fl36|) . and the parallel component of the fluid equation, Eq. fHTl) . respectively 

(MLg) - Ap) {~Kab + QabK) + Aii'(p + p)gAB - >^9ab'^nP = -]:(^9ab , (55) 



and 



{p + p)\/aK + K\/aP = . 



(56) 



In this case, the extrinsic curvature is intertwined with both the energy density and the 
pressure of the fluid via the non-minimal coupling A, being in addition sourced by the brane 
tension cr. As in the case described next, the role of the brane tension seems superfluous 
for generating the discontinuity of the bulk geometry at the position of the brane and thus 
accounting for the singular presence of the brane in the bulk space. Although the system 
is intricately entangled, we have just enough equations and constraints to be able to solve 
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it unambiguously given initial conditions. The procedure follows that suggested above for 
the general system. For further insight into a possible nature of such bulk field, we draw 
the attention to the next case. 



B. Cosmological Constant Fluid: p + p = 



Another case contemplated by Eq. (H2|) is that of the fluid equation of state induced on 
the brane being p + p = 0, which corresponds to the bulk fluid inducing a cosmological 
constant on the brane. The {AB) and {NN) components of the induced Einstein equations 
become respectively 



Gab - KKab + i^Qab [KcdKcd + K^) 



+ 



(^P(5) - Ap) 



A {uaub + Qab) KVn{p + p) 
-9ab + A (VaVb - gAB^^h) 



+ gAB^(5) 



p = 



(57) 



{M' - Ap) - 



-R + K'{--1] -KcdKcd 



+ A 



(5) 



+ XulK^Nip + p) + 



p-^Ka=0 



(5^ 



whereas the parallel component of the fluid equation and the stress-energy conservation 
condition become respectively 

{1 + X{R- 2K^ - KcdK'"'' + 2KVn}) [{uaUb + 9ab) Va(p + p)] 
- XK'^ {uaUb + Qab) ^a{p + p) 



+ A 



2 {uaUb + Qab) ^aK + 2KVa {uaUb) " M ^ + ^ ) '^bUn 



Viv(p + p) 



- [2XK {Kab - QabK) + qab {1 + X{R-K^~ KcdK"""])] VaP = 



(59) 



^{d - 1) (M|,(5) - Ap) VbK - X (^Kab - \9abK^ VaP 
XK {uaUb - UcQab) Va(p + p) = . 



The corresponding matching conditions are 



(Mp(5) - Ap) {-Kab + QabK) - XqabVnP = -t^cfQab 



(60) 



(61) 



and 



VaP + uaubVb{p + p) = 
14 



(62) 



We now proceed to investigate this case in more detail. Since p + p = everywhere on 
the brane, then from boost and translation invariance we must have that Va(p + p) = 
everywhere on the brane. Since both p+p and its first derivative along the parallel directions 
to the brane vanish on the brane, then so must vanish its second derivative. However, terms 
in Vn{p + p) or VArVA(p + p) do not necessarily vanish. From Eq. (162|) it follows that 
VaP = ^ aP = 0, and from Eq. ( !60l) that VbK = 0. Moreover, since V^(p + p) = 
implies that V^p = — V^p, whereas Vp = Vp implies that V^p = V^p, we must have that 
V^p = V^p = 0. We can then solve the Einstein, Eq. (1571) . and the fluid equation, Eq. ( 159|) . 
iteratively for the evolution of the induced metric and Vat (p+p) on the brane, with the sole 
ambiguity residing in the fluid velocity in the bulk. The value of the extrinsic curvature 
will then be given by Eq. fl58|) . We also note that with the bulk field behaving on the brane 
as an effective cosmological constant the role of the brane tension becomes superfluous. 

Furthermore, should we allow for Va{p +p) 7^ 0, then the discontinuity of the extrinsic 
curvature across the brane would be further sourced by the evolution of p + p on the brane 
and thus generated dynamically according to the equation of motion for the fluid induced 
on the brane. This would be an interesting idea to pursue further. Despite the increased 
complexity of the problem, the coupled system would still be solvable. 

C. Vanishing Extrinsic Curvature: K = 

The remaining case is that of a vanishing extrinsic curvature, where the {AB) and [NN) 
components of the effective Einstein equations become 




and the effective equations for the fluid 



Va [(p + p) {uaUb + gAs)] - Qab^aP = 



(65) 



(p + p) {uaUb + Qab) ^aK = , 



(66) 
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where we used the Israel matching condition 

AV^P = \(y . (67) 

For K = 0, the brane tension is supported by the discontinuity in the energy density of 
the bulk fluid across the two sides about the brane. The matching condition for the fluid 
equation then yields 

{p + p)ubUn {I + \R) = Q (68) 
and consequently R = —1/X.^ Then Eq. flM|) reduces to 

V^p-^(^^ + A(5)j = 0. (69) 

This equation can be solved for p given initial conditions at the intersection of the brane 
with the bulk past inflnity, obtaining the evolution on the brane of the energy density of 
the bulk perfect fluid in terms of the parameters of the bulk space. The solution must also 
be subject to the reproduction of a consistent bulk cosmological constant, in the case of 
anti-de Sitter A(5) < 0. Upon substitution of Eq. flM|l . Eq. (I63l) becomes 

(M|,(5) - Ap) {gab - ^9Ai^ + AV^VbP = . (70) 

From Eq. (l66i) . we must have that V aK = 0. We can then use the solution for p to solve 
Eq. fITOl) for qab and Eq. fl65l) for p given ua- Hence, ioi K = the system decouples 
and can be solved straightforwardly. This case is, however, too limited since it does not 
constrain the evolution of Vivp capable of generating the brane tension a in Eq. fl^TI) . 



V. CONCLUSIONS 



In this paper we have considered a modifled gravity model where the curvature scalar 
couples non-minimally to the matter Lagrangian density, which here we realize for the case 



of a perfect fluid. As discussed in Ref. j4|, in four dimensions this model can potentially 
account for the observed rotation curves of galaxies without recourse to dark matter and 



^ The other possibiUty, p+p = 0, would be but a particular case of the bulk fluid behaving as a cosmological 
constant on the brane, as discussed previously. 
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suggests a solution to the Pioneer anomaly. In addition to this coupling, we have also con- 
sidered an extra spatial dimension, so that our spacetime is embedded in a five-dimensional 
bulk space where gravity is described by a five-dimensional Einstein equation. 

We find that the resulting model is well defined for the considered physical variables 
and can be solved for given initial conditions. The new terms that arise from the bulk- 
brane decomposition yield quite interesting consequences. We found three particular cases 
which conform to the matching conditions upon the assumption of Z2-symmetry about 
the position of the brane and investigated their contribution to the intricately entangled 
system of equations of motion therein induced. These cases are the following: a) refiecting 
boundary conditions on the brane, ua = VatMtv = 0; b) an induced fiuid equation of state 
characteristic of a cosmological constant, p + p = 0; and c) a vanishing extrinsic curvature, 
K = 0. 

Both ) and b) seem to render the presence of a brane tension a superfluous for 

generating the discontinuity of the bulk geometry at the position of the brane, where the 
energy density p and the pressure p of the fluid can source the discontinuity of the extrinsic 
curvature in a dynamical manner. 

Furthermore, case b) can be regarded as a generalization of a cosmological constant 
scenario, where the bulk fluid induces on the brane a cosmological constant capable of 
supporting the presence of the brane. This implies that the evolution on the brane of 
the energy density of the bulk perfect fluid determines the behaviour of the cosmological 
constant term on the brane. It is well known that an evolution in terms of cosmic time 



as p oc t"^ is consistent with the value of the vacuum energy density at present [l6|, 113]. 
Furthermore, this positive contribution for the brane cosmological constant can have impli- 
cations for inflation and for the late time acceleration of the universe. This contribution can 
also have a bearing on the cosmological constant problem since the natural background for 
fundamental theories, such as supergravity and superstring/M-theory, is the anti-de Sitter 
space, which on the brane requires a compensating de Sitter contribution. The feasibility 
of a scenario along these lines will be considered elsewhere. 

Finally, case c) allows for the decoupling of the system of equations, with the energy 
density relating with the bulk cosmological constant A(5). The presence of the brane is 
supported by the interaction of the brane tension with the discontinuity across the brane 
of the fluid energy density, which on the brane is governed dynamically by a completely 
deflned equation for p. 
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APPENDIX A: TENSOR DECOMPOSITION IN GAUSSIAN NORMAL CO- 
ORDINATES 



Here we derive the projection of the Einstein tensor along the parallel and normal direc- 
;ions to the brane, namely Gab, Gan and Gat at, using the Gaussian normal prescription 

HI 

^'^Ve-^es = T^BA^c - KabCn , (Al) 
(^Ve-^CA = K^ec , (A2) 



(5)v- 

V p 



K^ec . 



(A3) 



The Gaussian coordinates split the bulk space into a four-dimensional spacetime times a 
one-dimensional space, so that 



/ 



(5), 



\ 



Here, A, B,... denote directions parallel to the brane, while N denotes the orthogonal 
direction. 

We aim to obtain decompositions of the form 

^^^Gab = terms without extrinsic curvature -|- terms with extrinsic curvature , (A5) 

so that the term without extrinsic curvature will be identified with ^^^Gab- We use the 
Einstein equations in the bulk to eliminate ^^^Gab and thus obtain the Einstein equation 
parallel to the brane. However, this equation will still contain terms in the extrinsic cur- 
vature, which can be related to other quantities on the brane through the Israel matching 
conditions. In this calculation we use geodesic coordinates on the brane, i.e., on the metric 
^^^qab- However, this does not mean that the connection components V^ab can be ignored, 
since their derivatives are non-vanishing. 
The Riemann tensor can be written as 



(A4) 



R{Y,W)Z = VyV 
where F, W and Z are vector fields, so that in components: 
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(A6) 



(A7) 



On its turn, the Ricci tensor is given by 



R^p = dx^iR^^^^e,) . (A8) 

Let us now compute the decomposition of Rab- First, one computes R'' acb(^j and 
R'^ANB^J- Noticing that the indices J and L refer to all five dimensions, it follows that 

^'^R'aCBCj = (^Ve-e ^'Ve-,e4 - {C ^ B) 

= (^Ve-e {T%eD - KabCn) -{C^B) 

= ^ab,c^d - KAB,ceN - KabK^bd -{C^B) . (A9) 

The symbol = denotes equality only for the case of geodesic coordinates on the brane. We 
also note that 

(^Ve-c {KABdx^ ® dx^) = dc (Kab) dx"^ ® dx^, (AlO) 
which is equivalent to Kab\c = Kab,c-, where \C denotes ^^^Vg-^. Then, 

^'^R'acbCj = {r^^c - ^ac,b) + {KacKE - KabK^) eo 

+ {Kac\b - Kab\c) ejv ■ (All) 

Computing now R'^anb^J ^^^^ 

^'^R^ANBej = ^'^Ve-, {r'^ABeD - KabSn) - ^'^^es ^'^^S^Sa 

= r%^ei, - KAB,NeN - K^iBeo + K^KdbBn ■ (A12) 
The five-dimensional Ricci tensor components can now be written as 

^'^Rab = dx"^ ( ^'^R'albCj) = dx"" ( '^'^R\cBej) + dx"" ( ^'^R'anbCj) 
^ ^AB,C ~ ^AC,B + KacKb — KKab — Kab,N + K^KcB 

= - KKab + 2X^0^^^ - Kab,n , (A13) 

where K — Kab ''^^g^^- Moreover, we compute 

(^Ve-^ {KABdx'^dx^) = KAB,Ndx^dx^ + KABdx^ ^^^Vg^dx^ + KABdx^ (^Ve-^(ix^ 

= KAB,Ndx^dx^ - 2KAcK^dx^dx^ , (A14) 

which yields 

Kab;N = Kab,n - 2KacK^ . (A15) 
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Finally, '^^^Rab can be rewritten as 

= '^^^Rab - KKab - Kab;N . (A16) 

Analogously, the other components of the Ricci tensor are given by 

^'^Rnn = -KcdK^''' - , (A17) 
= - K\A , (A18) 

so that the curvature scalar is expressed as 

(^)i? = g^"" '^'^Rab + ^''^Rnn = '^'^R -K'- KcdK^'' - 2Kn . (A19) 

Then, the decomposition of the five- dimensional Einstein tensor results as follows 

(5)^ _ (5)d 1 (5)^ (5) D 

<^AB — ^AB - 2 9AB 

= ^^^Gab + 2KacK^ - KKab - Kab,n 

+ h^)g^^{KcDK^'' + K' + 2K,^) , (A20) 

(^)G^^ = ^'^R^A - \ ^'^QNA ^'^R 

= K^j, - K\A , (A21) 

= -^i^)R+^-{K'-KcDK^^) , (A22) 

Equations flA20l - [A22p are known as the Gauss-Codacci relations. 

In order to derive the matching conditions in our model, we point out that the Einstein 
equation, Eq. ([7]), contains second derivatives of the energy density p. We need to express 
them in Gaussian coordinates. Under a transformation of coordinates ej = AjCj,, where 
A J are the components of the versor ej expressed in the e,^ basis, 

V;,v,p a^a:;V;,v,p 

= A^v^ (a:;v.p) - A^ (v^a:;) v.p 

= ViVjp-V^jjVlP . (A23) 

Then, terms of the form V^Vi^p decompose along the directions parallel and orthogonal to 
the brane as follows 

(AB) VaVbP + KabVnP , (A24) 

{AN) VaVnP - K^VbP , (A25) 

{NN) VnVnP . (A26) 
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Consequently, 

(5)np = np + VivVTvP + i^VjvP • (A27) 
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